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Based on the method of matched asymptotic expansion 
and on a time-dependent variational analysis, we study the 
dynamics of a vortex in the large-condensate (Thomas-Fermi) 
limit. Both methods as well as an analytical solution of the 
Bogoliubov equations show that a vortex in a trapped Bose- 
Einstein condensate has formally unstable normal mode(s) 
with positive normalization and negative frequency, corre- 
sponding to a precession of the vortex line around the center 
of the trap. In a rotating trap, the solution becomes sta- 
ble above an angular velocity f2 m characterizing the onset of 
metastability with respect to small transverse displacements 
of the vortex from the central axis. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40. Db 



Recent experimental observations of Bose-Einstein 
condensation in trapped alkali-atom gases at ultra-low 
temperatures have raised the question of vortex 

states similar to those in superfluid helium. Unlike su- 
perfluid helium, however, the harmonic trap makes the 
alkali Bose condensates nonuniform. 

Numerical studies of a vortex in small to medium ax- 
isymmctric condensates have determined the condensate 
density , along with the low- lying excitation spec- 
trum [|6|]7]] . These analyses found a nontrivial eigenmode 
with positive norm and negative frequency, and a varia- 
tional treatment || suggested that such a system might 
be unstable. In the Thomas-Fermi (TF) limit, the pres- 
ence of a vortex [PHTll generally splits the eigenfrequen- 
cies of a vortex-free condensate by a small amount of 
order (d/R) 2 <C 1 (but these analyses ignore the possibil- 
ity of the nonperturbative negative-frequency mode that 
occurs only when a vortex is present). Separate studies 
have considered vortex formation by a moving object ]T^ ] 
or a rotating force [^3| . 

The existence of a positive-norm eigenstate with 
negative-excitation energy |(^,|7j formally implies an in- 
stability of the vortex state because the condensate itself 
corresponds to zero energy. The present work obtains 
explicit analytical negative-energy solutions for a vor- 
tex in a large condensate, holding when the interatomic 
and trap interactions predominate and the kinetic en- 
ergy of the nonuniform density is small. This TF limit 
is valid for R/d » 15(iVa/d) 1 / 5 > 1 Q, where R is the 
mean condensate radius, d = (dxdydz) 1 ^ 3 is the mean 
harmonic oscillator length [here, d 2 — h/Mu)j, with M 
the atomic mass and u>j the corresponding trap frequency 
(j = x,y, z)], N is the number of atoms in the conden- 
sate, and a > is the s-wave scattering length. Typi- 



cally [jl5| d ~ a few /im, a ~ a few nm, and the exper- 
imental value N > 10 6 amply satisfies the TF criterion 
that N > 10 3 . In the TF limit (R > d), the vortex core 
radius £ ~ d 2 /R is small compared to both d and R. 

In this paper we assume that most of the particles 
are in the condensate, neglecting scattering with ther- 
mal background atoms, as confirmed by recent experi- 
mental studies of the lowest-lying normal modes (jHj 17 
for low temperatures. A condensate containing a g-fold 
quantized vortex in an axisymmetric trap is known to 
have only stable normal density modes ex e tm ^ for m = 
and \m\ > 2\q\ (so that unstable modes occur only for 
< |m| < 2|g|) jl(|. Here, we exhibit unstable solu- 
tions with \m\ = 1 that have positive norm and negative 
frequency in which the vortex line precesses around the 
center of the trap. 

Dynamics of a vortex in the trapped condensate — Let 
us consider a condensate in a nonaxisymmetric trap that 
rotates with an angular velocity Q around the z axis. At 
zero temperature in a frame rotating with the angular ve- 
locity £lz, the evolution of the condensate wave function 
'F is described by the time-dependent Gross-Pitaevskii 
equation 18 ,19]: 
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where V tT = \M (iv 2 x 2 + ui 2 y 2 + ^ 2 z 2 ) is the external 
trap potential, g — Airffa/M > is the effective inter- 
particle interaction strength, (J,(£l) is the chemical poten- 
tial in the rotating frame, and <f> is the azimuthal angle in 
cylindrical polar coordinates. We assume that the con- 
densate contains a g-fold quantized vortex parallel to the 
z axis located near the center of the trap at the trans- 
verse position p (t). In this section we use the method 
of matched asymptotic expansion to determine the vor- 
tex velocity as a function of the local gradient of the 
trap potential and angular velocity f2, generalizing two- 
dimensional results obtained by Rubinstein and Pismen 
p0| ] to the case of a three-dimensional rotating harmonic 
potential. The method applies when the external poten- 
tial does not change significantly on distances comparable 
with the core size \q\£ <C R± (TF limit); it matches the 
outer asymptotics of the solution of Eq. ([!]) in the vor- 
tex core region (\p — p \ < \q\0 with the short-distance 
behavior of the solution in the region far from the vortex 
core (\p — Pol ^ Recently the method was used to 

analyze the drift of an optical vortex soliton in a slowly 
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varying background field [|T) . 

To find the solution in the vortex core region, one 
may consider Eq. (Q) in a coordinate frame centered on 
the vortex line that moves with the vortex velocity V 
(Viz). The solution is assumed to be stationary in the 
comoving frame and satisfies the equation: 



2M 



V 2 + Vtr + s|*| 2 - a*(A) + itMdj 



-my ■ v i * 



(2) 



In the vortex core region we may seek a solution in the 
form of an expansion in the small parameter £/R±: 



# = [\V (p,z)\ - X (p,z) cos <f>]e 



iq<f>—ir}(p,z) sin ^ 



(3) 



where is the condensate wave function with Vtr re- 
placed by V z — \MlJ\z 2 and \ and rj characterize the 
perturbation in the absolute value and phase. Physically 
is the analogous wave function for an unbounded con- 
densate in the xy plane with the same chemical potential 
The polar angle <f> is measured from the direction 
of the local potential gradient Vj_Vt r (p )7 an d P is the 
local radial cylindrical coordinate. For p ^> \q\^ the per- 
turbations have the asymptotic form: 
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Here we omit terms containing f2 explicitly because they 
are proportional to the small parameter hfl/ p. The pa- 
rameter A must be determined by matching the solutions 
(^) with those far from the vortex core; in fact, A depends 
on f2. 

To lowest order in the small parameter £/-Rj_, Eq. (|l|) 
far from the vortex core reduces to an equation for the 
condensate phase only (* = |*|e ): 
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and p = p(Q) + hQq. 
Comparing the solution of Eq. (0) for p -C R± with the 
outer asymptotic form (0) of the core solution gives the 
value of the parameter A. With logarithmic accuracy we 
find 
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where is the mean transverse dimension of the con- 
densate. To find the velocity V of the vortex, one can 
use a solvability condition (Fredholm's alternative) of the 
first-order equation in the small parameter (see |20| for 
details). In the frame rotating with the trap, the vortex 
velocity is found to be orthogonal to the direction of the 
local gradient of the trap potential Vj_Vt r : 
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where p can be taken as the chemical potential for a 
nonrotating trap. The vortex follows an elliptic tra- 
jectory along the line Vtr = const. For q > and 
Q < fl m = 3\q\h (w 2 + w 2 ) ln(R_L/\q\£) /8p, the vortex 
moves counter-clockwise in the positive sense. With in- 
creasing rotation frequency fl of the trap, the vortex ve- 
locity (as seen in the rotation frame) decreases towards 
zero and vanishes at f2 = f2 m . At this rotation frequency, 
the effective potential for the vortex becomes flat near 
the trap center [see Eq. ( p8|) below]. For £1 > Cl m the 
apparent motion of the vortex becomes clockwise. 

From Eq. (]?]) it is straightforward to obtain the vortex 
position as a function of time: 



xq = eoRx sin (cot - 



yo = ±eo R y cos {wt + 0o ) 



where 



w = ± 
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This vortex motion near the trap center can be con- 
sidered a collective excitation of the condensate with 
an eigenfrequency u>. The two different frequencies 
given by ( |l0| ) correspond to solutions with different sign 
of normalization: in terms of the Bogoliubov ampli- 
tudes u(r),v(r), these solutions have the normalization 
/ d 3 r (|zi ± | 2 — |w ± | 2 ) = ±sgn(g). The solution with pos- 
itive norm has the energy 
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sgn(g) , 2 l 2 . 



3\q\h uj x U! y 
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If fl = 0, the excitation energy is negative and hence 
formally unstable. This negative energy implies the ex- 
istence of a state with lower energy that can become 
macroscopically occupied. Nevertheless, the condensate 
will be unstable only if there is a mechanism to transfer 
the system to the lower energy state |^2| , which requires 
the presence of a reservoir to conserve energy and angu- 
lar momentum. Thus a vortex in a nonrotating harmonic 
trap can remain stable at low enough temperatures when 
the noncondcnsate atoms are negligible. If the trap ro- 
tates, Eq. ([□]) shows that the solution becomes stable 
at |0| > Q m which coincides with the angular velocity 
at which the vortex at the center becomes metastablc 
(see below). Note that the precession frequency (10) in a 
nonrotating trap involves the product uj x LD y , whereas the 
metastable rotation frequency fl m instead involves the 
combination \{u x 
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Analysis of the Bogoliubov equations for an axisym- 
metric trap — The collective excitation energies E of the 
(in general, nonuniform) condensate are the eigenvalues 
of the Bogoliubov equations for the coupled amplitudes 
u(r), w(r) 



-|^v 2 + y tr + 2 3 |*| 2 - / ^))(; 



ihTLdtj) 



-g* 2 

-iMldct, 



(12) 



In this section we assume the trap potential to be ax- 
isymmetric, that is uj x — u> y = We can rewrite the 
Bogoliubov equations in the following form: 
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where 
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and 'J'o is the wave function for an unbounded condensate 
in the xy plane with the same chemical potential p(fl); 
its excitations obey the equation 
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For a g-fold quantized vortex at the center of the trap 
parallel to the z -axis, the unperturbed condensate wave 
function has the form * = e^^ol, and Eq. @ has 
the following exact pair of solutions [fill 



(-d p ± q/p) 
e-^ (d p ± q/p) , 



E Q = Thn, 



l*o|, (17) 
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where p is the radial cylindrical coordinate, and I 2 = 
4ir\q\ J dz |* (p = oo)| 2 « 16\/27r|g|^ 3 / 2 /35W 2 VM in the 
TF limit. These solutions are normalized so that 
/ d 3 r (|mq | 2 — \vq\ 2 ) — sgn They generalize those 

obtained by Pitaevskii |Q to the case of an axial trap 
with potential V z and describe the transverse motion of 
the vortex line as a whole relative to the condensate. It 
is obvious that the solutions U + and {/_ are orthogonal. 



We use Eqs. (|r^-(|l8|) as a zero-order approximation 
for a perturbation expansion of Eq. fll3| ) . The first-order 
shift of E requires the matrix elements (U±\V\U±). With 
logarithmic accuracy, their main contribution arises from 
the integration region \q\£ <C p <C R±, where \q\£ is the 
vortex core radius. For this region in the TF limit, we 
have 
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where 3|Wtf| 



\Muj 2 z z 2 and p = p(Q)+ hQq. The 
last two terms in ( |l9| ) reflect the gradient terms in the 
Gross-Pitaevskii equation and are small in the TF limit. 
Nevertheless, these small terms give the main contribu- 
tion to the matrix elements. The same expression ( pj| ) 
holds for the density of the unbounded condensate l^ol 2 
Thus the perturbation operator V 



but with lu± = 
becomes 



V=(-Mu; 2 ± {x 2 +y 2 ) + 
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and the off-diagonal matrix elements ^ U± \ V \ ) vanish. 
With logarithmic accuracy, the last term in Eq. (EC?) gives 



the main contribution to 



u±\v\u± 



2 , ,2 



d 3 r- 



(21) 



dp 



gP J p* 
gives a logarithmic factor- 



Integrating over 

ln(i?±/|g|0- 

As a result, we obtain two solutions for the excitation 
energy. They are normalized to the value sgn(=F<7), and 
the solution with positive norm has an energy 
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E = sgn(<7)/ifi 



3\q\h~ uj 
4M 



■ln(R±/\q\0- 



(22) 



where p can be taken as the chemical potential for a non- 
rotating trap. Formula (^2|) agrees with ([ll]) obtained on 
the basis of the method of matched asymptotic expansion 
(for uj x = uj y = ujj_). 

Energy of a vortex in a Bose-Einstein condensate in a 
rotating trap — In a frame rotating with angular velocity 
VLz, the energy functional of the system is 



dVl — 



ll*l 4 - 



(23) 



Consider a g-fold quantized vortex displaced from the 
center of the trap, with transverse coordinates x = xq, 
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y = Ho. In the TF limit, the vortex-induced change in 
the condensate density is negligible. Hence, the main 
contribution to the change AE(xo,yo) hi the energy of 
the system due to the vortex arises from the condensate's 
superfluid motion, and we can take |^| 2 in Eq. ( jgjj ) to 
be the density of the condensate without a vortex. A 
transverse shift of the Cartesian coordinate system places 
the origin on the vortex axis and yields 



|\|/| e «9<H-*So 



(24) 



9l*l 2 = ii - 2 M H( x + x o) 2 + rfiv + yo) 2 + ulz 2 \ , 

(25) 

where <f> is a polar angle around the vortex line and So is 
a periodic function of <f>. 

Varying the functional ( |23| ) gives an Euler-Lagrange 
equation for So . For a vortex-free condensate, we get 



So — — 



MVL ( IjjI - UJ 2 \ 

-j-i^-^l^ + ^iy + yo) ■ (26) 



If the condensate has a g-fold quantized vortex, then the 
contribution to the free energy from the phase factor iSo 
in Eq. (24) can be neglected with respect to the contri- 
bution from iq(f> for small angular velocity Q. Only for 
|Q| \q\K (ui 2 . + ujy) /4/i do the contributions from these 
two phase factors become comparable. For such large an- 
gular velocities, however, the function So is insensitive to 
the presence of the vortex. Thus Eq. (p3) with So given 



by Eq. (26) provides a good approximation for the con- 
densate wave function for general angular velocity fl. 

Substituting Eqs. (|24|)-(26) into the energy functional 
( p3| ) and integrating with logarithmic accuracy over the 
volume of condensate yield the following expression 
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(27) 



where R 2 = 2/j,/Muj 2 fixes the condensate's dimensions, 
R^ 2 = M(uj 2 x + u^)/4/z = \{Rx 2 + Ry 2 ) character- 
izes the mean transverse dimension, £ 2 = Ti 2 /2M(i, and 
no(0) = ix /g is the density at the center of the vortex- 
free condensate. For small displacements of the vortex 
(x <C R x , yo -C R y ), we have: 



&E(x ,yo) 
1 



fMiU 2 MO){g 2 ln(^ 



8q[iQ 



' ii 2 
) | Vo 

2 \^ "^^i; ^^"y 



5h {luI+uD 
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For |fi| > n c = 5\q\h(<4+ul)ia.(R ± /\q^)/8n = 
^\q\h/(MR 2 L )\n(R±/\q\^), the central vortex is sta- 
ble because AE(0, 0) is negative p3|. Furthermore, 
AE(xo,yo) has a local minimum at the origin for |fi| > 
Cl m = |^ci showing that the central vortex is metastable 
in the interval f2 m < |fi| < £l c . Since huj XtViZ <C \i in the 
TF limit, the thermodynamic critical angular velocity £! c 
for vortex formation and il m are both much smaller than 
the trap frequencies. 

Time- dependent variational analysis — Instead of solv- 
ing the time-dependent Gross-Pitaevskii equation (Q) di- 
rectly, one may consider a variational problem for the ac- 
tion obtained from the effective Lagrangian (we assume 
Q = in this section): 



L(t) 
2M 



d A r 



W| 2 - V t r(r)|*| 2 -±ff|*| 



(29) 



With a suitable trial function containing time-dependent 
parameters, the time evolution follows from the associ- 
ated Lagrange equations for these parameters. Although 
not exact, this method gives a clear physical picture of 
the system's behavior. It determined the low-energy ex- 
citations of a trapped zero-temperature condensate [p4[ , 
predicting the low-energy excitation spectrum of the con- 
densate and, in the TF limit, reproducing the formulas 
derived by Stringari p5| . In the present paper we use 
this method as an alternative way to study the normal 
modes of a condensate with a vortex. 

We are interested in the relative motion of a vortex 
inside the condensate, so it is natural to take a time- 
dependent trial function in the form: 

V(r,t)=B(t)f[r-p (t)]F[r-r, (t)} 

x Yl ex P [ixj<Xj(t) + iXjPjty)] . (30) 

j=x,y,z 

Here the function / (r) characterizes the vortex inside 
the trap and has the approximate form / (r) = e lq ^ for 
P ^ M£> the function F(r) describes the condensate 
density distribution and, in the TF limit, has the form 



F 2 (r) = 1 



y 



Rl{t) R*(t) Ri{ty 



(31) 



(28) 



The time-dependent vector T] {t) = (r)o x , rjo v , ?7oz) de- 
scribes the motion of the center of the condensate, while 
Po W = ( x 0i J/0) 0) describes the motion of the vortex line 
in the xy plane. The other variational parameters are the 
amplitude B, the width of the condensate Rj, and the set 
ctj, and (3j that describe the condensate's motion. These 
parameters are real functions of time, characterizing the 
macroscopic wave function of the condensate. Substitu- 
tion of the trial wave function pel) into Eq. (p9|) yields 
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an effective Lagrangian as a function of the variational 
parameters (and their first time derivatives). 

We focus on small oscillations, when the displacement 
Po(t) of the vortex from the trap center and the dis- 
placement T] (t) of the condensate are small relative to 
the width of the condensate; expanding for small p (t) 
and r) (t), we keep only the first nonzero corrections. 
The resulting equations for the motion of the vortex line 
(that is, for xq, t/o) ar e coupled only with r]o x ,r]oy and 
a Xl a y . This system of six first-order ordinary differen- 
tial equations can be solved explicitly. Apart from the 
well-known dipole-mode solutions that describe the rigid 
oscillation of the condensate and the vortex as a whole, 
these equations have another solution that corresponds 
to the motion of the vortex relative to the condensate. 
For this solution the vortex motion is described by for- 
mulas (||), (0), while the displacement of the condensate 
is given by 



VOy = T 



where 



2R n 



15e g 2 
2R X 



q\£j R x +R, 



■ sin (tot + (j> ) , (32) 



In 



R± 



R,, 



Rx + R, 



■ cos (cut + 4>o) , (33) 



3qhLU x ujy 

'- In 

4/i 



Rx 



(34) 



The quantity Xq/R x + y^/Ry = Eq remains constant 
as the vortex line follows an elliptic trajectory around 
the center of a trap along the line Vt r = const, and 
the energy of the system is conserved [as follows from 
Eq. (|27])]. These results agree with those obtained from 
the method of matched asymptotic expansion. The con- 
densate also precesses with relative phase shift tt at the 
same frequency, but the amplitude of the condensate mo- 
tion is smaller than that of the vortex line by a factor 
~ £ 2 In (Rx/\q\£,) /R x R y . As in the discussion of (10), 
these solutions have the normalization — sgn(w), and the 
positive normalization corresponds to the solution with 
negative frequency. In the present limit of logarithmic 
accuracy, the structure of the vortex core is inessential 
for the excitation spectrum. The precession frequency 
([m]) obtained in the variational approach coincides with 
formula (jl0|) (for ft = 0). 

In conclusion, we have used the method of matched 
asymptotic expansion (for a general three-dimensional 
rotating nonaxisymmetric trap), analysis of the Bogoli- 
ubov equations and a time-dependent variational pro- 
cedure to study particular normal modes of a trapped 
condensate containing a vortex. These modes have nega- 
tive energy and positive norm, with the vortex line mov- 
ing around the z axis. The excitation frequency in a 



rotating trap becomes positive when the angular veloc- 
ity f2 reaches the value that makes the central vortex 
metastable. 
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